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An analysis is presented of wave- vector dispersion in elliptically birefringent 
stratified magneto-optic media having one-dimensional periodicity. It is 
found that local normal- mode polarization-state differences between adjacent 
layers lead to mode coupling and impact the wave-vector dispersion and 
the character of the Bloch states of the system. This coupling produces 
extra terms in the dispersion relation not present in uniform circularly 
birefringent magneto-optic stratified media. Normal mode coupling lifts the 
degeneracy at frequency band cross-over points under certain conditions and 
induces a magnetization-dependent optical band gap. This study examines 
the conditions for band gap formation in the system. It shows that such a 
frequency-split can be characterized by a simple coupling parameter that 
depends on the relation between polarization states of local normal modes 
in adjacent layers. The character of the Bloch states and conditions for 
maximizing the strength of the band splitting in these systems are analyzed, 
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1. Introduction 

A number of publications have discussed the propagation of electromagnetic waves in one- 
dimensional birefringent stratified media. [1,2] Of particular interest is the formulation pi- 
oneered by P. Yeh that employs a translation matrix approach to discuss periodic linearly 
birefringent structures. The propagation of light across a single period is analyzed in con- 
junction with Floquet's theorem to determine the dispersion relation and Bloch waves of 
the system. This approach has been used to study the properties of layered media with 
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misaligned optic axes from one layer to the next. [2] Numerical solutions have been found 
for normal incidence of an optical beam for identical birefringent plates with alternating 
azimuth angles and layer thicknesses. [1] 

More recently the present authors have used this technique to study elliptically birefringent 
non-reciprocal media, [3] where elliptically polarized normal modes characterize the system 
locally. Analytic solutions were found for light incident perpendicularly into stratified media 
consisting of alternating magneto-optic (MO) layers with different gyration vectors [4] and 
birefringence levels. [3] Adjacent layers were assumed to have their anisotropy axes aligned 
to each other. This model captures some important features of one-dimensional magnetopho- 
tonic crystal waveguides, particularly the presence of magneto-optical gyrotropy and locally 
alternating birefringence levels. Such systems are presently being studied for use in integrated 
fast optical switches and ultra-small optical isolators. [5-7] Work on these systems extends 
prior theoretical and experimental efforts on magnetophotonic crystals in order to encompass 
the elliptical birefringence that often characterizes planar magnetophotonic structures. [8-12] 

A particularly interesting feature of periodic elliptically birefringent gyrotropic media con- 
cerns the character of the dispersion branch solutions to the Floquet theorem. The work of 
Pochi Yeh cited before [1] discusses the formation of a band gap away from the Brillouin zone 
boundary in linearly birefringent media with misaligned anisotropy axes. Merzlikin et.al. have 
pointed out that magnetically tunable band gaps can arise in stratified media that combine 
circularly and linearly birefringent layers. [13] In the present work we note that magnetically 
tunable band gaps can also exist in elliptically birefringent stratified magneto-optic media. 
A range of frequency splits is possible since local normal modes span a wide spectrum of po- 
larization states. The work presented here traces the origin of these band structural features 
to normal mode coupling arising from the simultaneous presence of gyrotropy and linear 
birefringence in the periodic system, as the Bloch states of the system propagate across layer 
boundaries. The dispersion relation for this type of system is found to contain additional 
terms that describe the formation of magnetically tunable band gaps away from the Brillouin 
zone boundary and that determine the magnitude of the gap. A new kind of parameter is 
identified herewith that characterizes the coupling of these layer-dependent normal modes. 
The underlying phenomenon is the continuity of transverse electric and magnetic field com- 
ponents across inter-layer boundaries. Such a coupling can yield magnetization-dependent 
optical band gaps inside the Brillouin zone as discussed below. The formulation presented 
herein can serve a tool for the design of magnetically tunable band gaps in non-reciprocal 
magnetophotonic structures. 

After introducing the formalism to be employed in this work (Sec. |2]), it is shown that the 
transfer matrix for elliptically-birefringent magnetophotonic crystals can be parametrized 
in terms of an inter-modal coupling parameter (Sec. [3]). The conditions for the frequency 
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splitting of degenerate Bloch states is then discussed (Sec. H]) and the dispersion relation 
for magneto-optic layered structures in the presence of elliptical birefringence is derived and 
discussed (Sec. [5]). Section [6] examines the character of Bloch states for these systems. 

2. Waves in a birefringent magnetophotonic medium 

In the optical wavelength regime, the permeability of a birefringent uniaxial magnetooptic 
medium is very close to the permeability of vacuum /iq, its relative permeability close to 
unity. The relative permittivity tensor e of the medium for magnetization along the z— axis, 
has the form 

6 ^^xy ^yy 

V e,, 

where we assume no absorption of the light in the medium. This implies that all components 
of the relative permittivity {eij,i,j = x,y,z) are real and it is not assumed that e^^^^ = eyy. 
By solving the wave equation upon normal incidence of a monochromatic plane wave (with 
time dependence exp (iut)) propagating parallel to the z axis on a birefringent magnetooptic 
medium, one obtains eigenmodes 

^ cos a ± sin a ^ 
±i cos a — i sin a 

\ / 



(2) 



corresponding to the refractive indices 



ni = -e±^A^ + ely. (3) 

Here e = {eyy + exx)/2, A = {eyy — exx)/2, and tan(2Q!) = A/e^y The propagation constant 
in the medium is defined as 

(3± = -n± (4) 
c 

where c is the speed of light in vacuum. 

3. One-dimensional birefringent magnetophotonic crystals 

Consider a plane wave normally incident on a periodic stack structure of alternating ellipti- 
cally birefringent magnetophotonic layers (Fig. [1]), where the elliptical birefringence param- 
eters of adjacent layers may differ, but the anisotropy axis are aligned. The Bloch states for 
this system can be expressed in terms of local normal modes. Thus in the n-th layer the 
optical electric field can be written as 

E(")(z) = (E(^)e^'^+'(^-^") + E£)e-^^+'(^"^"))eJ) (5) 
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Here E, 
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, 4) are the complex amplitudes of the partial waves corresponding 



to each normal mode. The Bloch states for this system satisfy the Floquet-Bloch theorem 
through the following eigenvalue equation. 



(6) 



where the transfer matrix '^(^-'^,^+1) Yelates the four eigenmode amplitudes E in the second 
layer of a unit cell (the region between z = {n — 2)A and z = [n— 1)A) to the corresponding 
amplitudes in the second layer of the adjacent unit cell (the region between z = {n — 1)A 
and z = nA). K is the Bloch wave vector and A is the period of the periodic structure. With 
this functional basis, the transfer matrix 'j^in~i,n+i) acq^j^gg following form: 
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Here d^'^^ is the thickness of layer n in the stack and = — The other 

elements of the T*^""^'""*"^) matrix can be obtained by the following symmetry relations 
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The 

4. Normal mode coupling and the splitting of degenerate Bloch states 

Pochi Yeh has discussed the propagation of electromagnetic waves in alternating linearly 
birefringent layers where the normal modes differ in adjacent layers due to anisotropy axes 
misalignment. [1] He points out that in this case a new type of constructive interference in the 
scattered waves arises from the coupling between slow and fast waves. Forbidden frequency 
zones or band gaps can now appear away from the Brillouin zone boundaries. He calls this 
an exchange Bragg condition because forward propagating fast (slow) Bloch states couple to 
backwards propagating slow (fast) Bloch states. A similar type of constructive interference 
exists in elliptically birefringent magneto-optic layered systems with different local modes in 
adjacent layers. It is shown here that such local normal mode variations lead to the opening 
up of a band gap in the Brillouin zone. This effect can be traced to the presence of normal 
mode coupling between adjacent layers. 



Notice that the unit cell transformation matrix x^""^'""'"^^ Eq. ([7]), depends on the relative 
elliptical birefringence parameter and the individual normal mode propagation con- 

stants for each layer. When the parameter equals zero the normal modes in adjacent 

layers are the same and they remain uncoupled. This can be seen explicitly from the form 
of the X^""^'""*"^) matrix formulated in terms of normal modes. It is clear from the expres- 
sion that if = the off-block-diagonal components of the x^""^'""*"^) matrix are zero 
and there is no admixture of the local normal modes. This situation is similar to the case 
of a periodic layered medium consisting of isotropic layers, where transverse electric (TE) 
and transverse magnetic (TM) waves, or right- and left-circularly polarized waves remain 
uncoupled. [10] 

When differs from zero the Bloch states require an admixture of local normal 

modes, since the off-block-diagonal components of x*^""^'""*"^^ differ from zero. Moreover, the 
strength of the coupling, parametrized as the weight of the off-block diagonal terms, can 
be seen to increase as sin 2x''"'"'"'"^\ We thus see that parametrizes the degree of 

admixture of the normal modes. In the next section we shall discuss the form taken by this 
coupling in the wave-vector dispersion of the system. 

As chan ges away from zero the polarization state of the normal modes changes 

according to Eq. ([2]). Local normal modes in adjacent layers acquire different polarization 
states. These normal modes in adjacent layers are now coupled by the continuity of tangential 
components of the magnetic and electric fields across the boundary. Changes in normal mode 
polarization thus affect the solution to the Floquet-Bloch theorem through mode coupling 
across the boundary. 

This effect can be seen in Fig. (j2]) below. The new band gap that develops away from the 
Brillouin zone edge will be denoted as Yeh band gap since it was first noted by P. Yeh for 
the case of linearly birefringent stacks. Although the case we are considering differs from 
Yeh's treatment due to the elliptical birefringence and non-reciprocity of the magneto-optic 
systems under consideration, it is still the presence of birefringence that is the key to the 
formation of this new type of gap. 

5. Dispersion relation as a function of inter-modal coupling parameter 

In this section we explicitly find the dependence of the dispersion relation on the inter-modal 
coupling parameter The derivation and final form of the dispersion relation highlight 

the dependence. We will consider the polarization state dependence on in 

the next section. 

If one maintains the same index contrast in the photonic crystal and the same normal 
mode propagation constants in each layer but allows the polarization state of the normal 
modes to vary, the only parameter that changes in the X^""^'""*"^-' matrix is and the 
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band splitting that emerges is a result of normal mode coupling. 

The characteristic equation of the T'^""^'""'"^^ matrix in Eq. ([71) is a polynomial function of 
A of order four. A denotes the eigenvalue exp (iKA) of the Floquet-Bloch expression Eq. ([6]). 
It does not denote wavelength as is otherwise customary. For MO materials the s are 
much smaller than fi/s in the T*^""^'""^^^ matrix. The characteristic equation can be written 
as: 

^(A) = /(A) + h{X) sin^ + j(A) sin^ 2x^"'''+^^ 

+k{X) sin^ 3x("'"+^^ + /(A) sin^ 4x^"'"+^^ = 0, (15) 

where /(A) is the characteristic equation of T*^""^'""*"^-' matrix when = 0, h, j, k, and 

/ are functions of n^J!l''^~^^\ j3^'^~^^\ and We note that for an MO periodic structure, 

f^h^j^k^l. The zeros of /(A) are denoted by Aq (which corresponds to the case 
where the off-block diagonal elements of T*^*^"^'""*"^) equal zero) and of g{\) by Aq. Upon first 
order expansion of g{X) around Aq we obtain 

g{X) = g{X'^)+g'{X'^){X-X',) 

^ h{X',) sin^ + ifiK) + h'iX'o) sin^ x("'"+^))(A - A[,). (16) 

From this expression one can obtain Aq in terms of Ag 

h{X'Q) sin2 



(Ao - K) 



/'(A^) + /i'(A'o)sin2x("'"+i) 



Aq ~ ° ~ /'(A') ^ ' y^') 

where a prime on the functions means the first derivative with respect to their argument. 
In Eq. ( ITT]) use has made of the fact that 5'(Ao) = and /(Ag) = 0. Terms proportional 
to j, k, and / and their derivatives have been neglected as these are all small polynomial 
expressions in A. From this simple analysis we can see that the difference in eigenvalues is 
directly proportional to sin^ x''"'"'''^^- 

Let us now consider the dispersion relation Eq. (35) in Ref. 3 for the case where ;^("'"+^) = 
0, given by 



cosi^iA = ( cos(/3l"+')rf("+^)) cos(/5i")d(")) - ^iV± sin(/5i"+')ci("+i)) sin (/?£") d^"); 



(1^ 



In this case band gaps appear only at the boundary of the Brillouin zone, displaying complex 
solutions for KA. On the other hand, when differs from zero, the dispersion relation 
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acquires an additional term in the form of 'Si[h{\Q) / f (X'q)] sin^ ■^in,n+i)^ follows: 
cos K±A = 



cos(/3i")d(")) - ^iV±sin(/?i"+')d("+i)) sin(/3i"^d("): 



-3? 



0±. 



in) 



(19) 



For a 



Here 3? denotes the real part of its argument and N± = 
general MO material N± ^ 2. This extra term in the dispersion relation is responsible for 
band gap formation away from the Brillouin zone edges. A complex solution for Kj^A (the 
same treatment can be applied for K_A), and hence the existence of a band gap, occurs 
under the following conditions 

4 + 4u{X'qJ sin^ X^""'""^^^ 



cos 2/5A < 
cos 2/5A > 



-4 



'0+^ 

2 + A^+ 

4m(A[) ) sin^ 



(20) 



where pA = cc./(2c)(ni"^rfW + ni"+'^rf("+i)) and m(A[)^) = 3?(/i(A[,^)//'(A[,^)). 

If one maintains the same refractive index contrast and normal mode propagation con- 
stants for each layer of the periodic structure and allows the normal mode polarizations to 
change, the Yeh band gap width is a function of ;^("-'"+^) only. This band width increases 
monotonically with < ^("'"+^) < 7r/2, as the range of f3A expands according to Eq. fl20|) . 
From Eq. f|T9l) K acquires a complex solution only for negative u{X'q^) for a typical MO 
periodic structure. The upper bound of f3A will thus have larger values and the lower bound 
lower values as increases. Notice that these bounds occur under equality in Eq fl2U]) . 

This results in a wider Yeh band gap. 

As an example let us consider a model system composed of bismuth iron garnet (BilG) 
with typical values for the dielectric tensor in the near infrared region 

/ 6.5411 zO.018 \ / 



-in) 
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zO.018 
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and e("+^) 
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5.9859 zO.018 
-i0.018 6.1699 
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(-ZZ ) 



(21) 



This structure simulates the effective index variation of a magnetophotonic crystal on ridge 
waveguide. [14] In Fig. [21 we show the band structure for the case where x^"^'"^^^^ = in 
dashed lines, corresponding to the solutions for /(A). There is no band splitting in the band 
structure at the crossover point. On the other hand, when ^^('^■"+1) ^ 0, a band gap opens 
up, as shown by the solid lines. 

In Fig. [3] we show the variation of the Yeh gap bandwidth for the same structure upon 
variation of the off-diagonal components in the dielectric tensors, corresponding to the tuning 
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of magnetization by an external magnetic field applied to the photonic crystal structure. We 
also show in the same figure the variation of the Yeh gap bandwidth upon the variation of 
birefringence in the waveguide periodic structure. We maintain the same refractive index 
contrast in the periodic structure while the birefringence of the layers varies. In both cases 
the Yeh gap bandwidth increases with ;)(^{"'"+^). 



6. Bloch states in periodic birefringent media 

Let us define the n-th unit cell as the combination of layers of n and (n + 1). The translation 
matrix from the second layer of the (n — l)-th unit cell to the second layer of n-th unit cell 
is given by the x*^""^'""*"^) matrix in Eq. ([7]). Upon solving the eigenvalue equation ([6]), the 
eigenvectors are given by 



E 



K 



Ck 



Bn 
Cn 



(22) 



K 



where ck is an arbitrary constant. An, Bn, Cn, and Dn are functions of n^J^'^'^^^ , (3^''^^^\ 
and According to Eq. the Bloch wave solution in the second layer of the n-th 

unit cell is then given (up to a constant factor) by 



Er'^(^) 



-nA) 



giXAg(n+l) 



+ 

jKA-(n+l) 



(23) 



Note that this expression depends on the relative birefringence parameter through 
the components of the matrix x("~^'"+^). Figure H] shows the polarization states of the Bloch 
eigenmode at the interface between layer (n — 1) and layer n and at the interface between 
layer n and layer {n + 1) for artificially large values of the off-diagonal component of the 
dielectric tensors, in order to highlight the polarization variation of the Bloch state. In this 
example, we have taken the material dielectric tensors as 

/ 4 iO.l \ / 6 



~{n) 
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-zO.l 




iO.l 
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and e("+^) 



V 



-i2 




i2 
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\ 



ezz J 



(24) 



In general, the polarization state of the Bloch mode will also evolve as the wave transverses 
the unit cell, although for typical values of the dielectric tensors these changes will be small. 
Thus the character of the Bloch mode is affected by the coupling. 

Whereas Bloch states for non-elliptically birefringent gyrotropic one-dimensional stacks 
are still circularly polarized, elliptically birefringent stacks can have Bloch modes whose 
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polarization states differ from that of the normal modes in each layer, and that depend on 
^(n,n+i) according to Eq. fl23l). 

7. Conclusions 

Mode coupling as a result of periodic variations in the polarization state of local normal 
modes in elliptically birefringent non-reciprocal periodic structures is reported and discussed. 
Interlayer normal mode coupling in such media affects the polarization state of the Bloch 
waves and the wave- vector frequency dispersion. This interlayer coupling is absent in isotropic 
and uniform circularly-birefringent periodic media. As a consequence of local normal mode 
coupling, extra terms appear in the dispersion relation characterizing the formation of a 
frequency band gap inside the Brillouin zone away from the zone boundary. The band width 
of this gap is found to be parametrized by a single characteristic coupling constant, and 
is shown to increase monotonically with this coupling parameter. An expression for the 
latter is presented and shown to depend on the difference in diagonal components of the 
dielectric tensors and the gyrotropies of adjacent layers in the periodic structure. Thus a 
ready made tool for designing wavelength dependent band gaps in non-reciprocal periodic 
magnetophotonic structures and calibrating their band widths is presented. Bloch mode 
polarization states are found to differ from those of the local normal modes and to evolve 
into different elliptical states as the wave propagates down the crystal. These Bloch mode 
polarization states are found to depend on the strength of the coupling between local normal 
modes in different layers. 
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Fig. 1. Schematic diagram of a one- dimensional birefringent magnetophotonic 
crystal with period of A. The magnetophotonic crystal extend indefinitely 
in the x and y directions. A plane wave is incident normally to the layered 
structure. jalaliFl.eps 
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Fig. 2. The energy band diagram for BilG periodic structure (fifth branch) 
with (i^"-' = 0.4 and = 0.6. The dashed hne corresponds to the case 

where = in the transfer matrix. The sohd hne corresponds to a 

reahstic case with x = 0.14. jalahF2.eps 
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Fig. 3. The width of the Yeh band gap versus The Yeh band gap was 

calculated for the fifth branch of the band structure of the periodic structure of 
BilG with = 0.4 and = 0.6. In one case the average in the refractive 

indices is kept constant while is allowed to change (solid line). In the 

other, the diagonal elements of the dielectric tensors of adjacent layers are kept 
constant while the off-diagonal elements are allowed to change simultaneously 
(dashed line). jalaliF3.eps 
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Fig. 4. Polarization of Bloch wave traveling through a magnetooptic birefrin- 
gent periodic structure. The Bloch wave polarization is depicted on the bound- 
ary of each layer in a unit cell (solid ellipses) just before the Yeh band gap in 
the band structure of the medium. Dashed ellipses show the local eigenmodes 
polarizations e_|_ for each layer. jalaliF4.eps 
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